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ABSTRACT 

We investigate the impact of the intracluster medium on the weak-lensing shear power spectrum 
(PS). Using a halo model we find that, compared to the dark matter only case, baryonic pressure leads 
to a suppression of the shear PS on the order of a few percent or more for I > 1000. Cooling/cooled 
baryons and the intergalactic medium can further alter the shear PS. Therefore, the interpretation of 
future precision weak lensing data at high multipoles must take into account the effects of baryons. 
Subject headings: cosmology: theory — dark matter — galaxies:clusters:general — gravitational lens- 
ing — large-scale structure of universe 



1. INTRODUCTION 

Weak gravitational lensing is a promising technique for 
precision cosmology. Its chief advantage over other tech- 
niques comes from its sensitivity to the total density field, 
rather than just the baryonic component. Since the to- 
tal density is mostly due to the dark matter, which only 
interacts gravitationally, the statistical properties of the 
total density field can be calculated ab initio with much 
higher accuracy than can observables that are more di- 
rectly dependent on baryons. Combined with the enor- 
mous statistical power of proposed survey projects, such 
as the Large Synoptic Survey Telescope^ and Super- 
Nova/ Acceleration Probe^ , weak lensing is thought to be 
capable of determining (some) cosmological parameters 
to a few percent level (e.g. Refregier et al. 2004). 

Many studies of weak-lensing statistics assume baryons 
to trace dark matter exactly (e.g. Bacon, Refregier & 
Elhs 2000; Cooray, Hu & Miralda-Escude 2000; Hu & 
White 2001; Song & Knox 2003). This allows one to 
predict the weak-lensing shear power spectrum (PS) from 
the nonlinear dark matter PS that is well calibrated using 
TV-body simulations (Peacock & Dodds 1996; Ma & Fry 
2000; Smith et al. 2003). 

However, more than 10% of the total mass is in 
baryons, which do not follow dark matter exactly on 
small scales. Therefore the total mass distribution and 
its PS will deviate somewhat from those in a universe 
that has the same initial conditions except that all the 
mass is in dark matter.'^ Such a deviation in the mass PS 
will result in a difference in weak lensing statistics. Even- 
tually, it will be necessary to account for the baryonic 
influence on the shear PS as we reach for the projected 
statistical power of weak lensing surveys. 

Baryons exist mainly in three categories: (1) cool- 
ing/cooled baryons such as in stars and galaxies, (2) 
hot baryons or the intracluster medium (ICM), and 
(3) the intergalactic medium (IGM). The cooling/cooled 
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^ One may think of this as an A'^-body simulation of dark mat- 
ter that has the same parameters and initial mass PS (including 
baryonic features such as due to acoustic oscillations) as a full hy- 
drodynamical simulation. 



baryons alter the statistics of gravitational lensing 
(White 2004), because they can condense to a much 
denser state than dark matter, which in turn modifies 
the dark matter distribution (Kazantzidis et al. 2004). 
The ICM is smoother than cooling/cooled baryons, so 
that the effect might be much less pronounced (White 
2004). However, the ICM greatly out- weighs the baryons 
in the galaxies. Therefore, it could still have a significant 
impact on the shear PS. Indeed, we find at I < 3000 the 
ICM has a larger impact on the shear PS than do the 
cooling/cooled baryons. 

For completeness we note here that the IGM is much 
more diffuse than the other two, but it is the largest 
reservoir of baryons. Since the density fluctuations on 
the angular scales we consider are dominated by high- 
density and massive objects, we neglect the IGM for the 
present, but plan to study its effect using hydrodynami- 
cal simulations in the future. 

In this Letter we demonstrate the effect of hot baryons 
on the weak-lensing shear PS. Assuming an isothermal (3- 
model (Cavaliere & Fusco-Femiano 1976), we are able to 
extend the halo model (for a review, see Cooray & Sheth 
2002) to include hot baryons in the total mass PS, which 
is then used to calculate the shear PS. Aside from the 
uncertainties of shear statistics due to halo parameters 
(Takada & Jain 2003), this method is not sufficiently ac- 
curate for interpreting future data, because we have not 
included all the baryons, and because we have made sim- 
ple associations of the ICM profile with its dark matter 
halo. Nonetheless, it provides a framework for analyti- 
cally estimating the effect of uncertainty in the clustering 
properties of hot baryons. 

2. MASS POWER SPECTRUM 

Realistic halo profiles (e.g. Navarro, Frenk & White 
1996; Moore et al. 1999) enable the halo model to repro- 
duce the nonlinear dark matter PS. With the one-halo 
PS Pdm(fc) and two-halo PS Pf^{k), the dark matter PS 
Pdm{k) can be approximated by (Peacock & Smith 2000; 
Ma & Fry 2000; Seljak 2000; Scoccimarro et al. 2001) 

Pdm(fc)=Pdm(fc)+Fdm(fc), (1) 
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where pd is the mean density of the universe at the 
present, Ud[k , md{i')] is the normahzed profile of a dark 
matter halo with mass rrid in Fourier space (e.g. Cooray 
& Sheth 2002), and Pl is the linear mass PS. The dark- 
matter density profile is cut off at the virial radius. The 
height of the density peak, i/, is defined as 



V = [Sc{z)/a{md)Y 



(2) 



where S^{z) ~ 1 .686n^°°^^ {z) (Hiotelis 2003) is the over- 
density of a spherical region that collapses at redshift z, 
fljj^^z) is the ratio of matter density to critical density 
at z, and CT(md) is the rms value of the density contrast 
5pd/pd within a radius of (3TOd/47rpd)^^'^ at z. Equation 
(2) also defines md(z^). The functions /{v) and &(z^), re- 
lated to the halo mass function and bias, respectively, 
have the forms (Sheth & Tormen 1999) 

uf{u)=A{l + u-P)uy^e-''^/^ 

b{i^) = 1 + (i^i - 1) + 2pd-\l + <)-\ 

with i/i = 0.707;^ and p = 0.3. The normalization con- 
stant A satisfies the constraint / f{i')di'= 1. 
Equation (1) can be extended to include hot baryons: 



P{k)=J2fifiPijik), 



(3) 
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= j /(z^)dz/ md{i')ui[k,md{i')]uj[k,md{i')]/ pd, 
P^ = Pi.{k) J f{u)b{iy)diyudk,m^{i.)] 

X J f{u')b{iy')du'uj[k,md{u% 

where the subscripts i,j = d for dark matter and b for 
baryons. The coefficients, fd and /b, are the fractions of 
dark matter and baryons in total halo mass, respectively. 
Obviously, Pij{k) = Pji{k), fd + fb = 1, and equation 
(1) is a special case of equation (3) when /b = 0. 

For equation (3) wc have assumed a fixed gas-mass 
fraction /b. In principle, there is a distribution of /b for 
a given halo mass, and this distribution can be dependent 
on halo mass (or on the ICM temperature, see Arnaud & 
Evrard 1999; Mohr, Mathiesen k Evrard 1999, hereafter 
MME99). Although the dependence is found to be weak 
for clusters (Ota & Mitsuda 2004, hereafter OM04), it 
may not be true for low mass systems. Hence, an average 
over the distribution should be incorporated to calculate 
the PS. We assume for convenience that /b equals the 
cosmic mean baryon fraction of 0.13. 

The Fourier-space baryon profile Ub(A:,md) is calcu- 
lated from the /3-model, pb{r) = /9bo(l + r^/r^)"^^/^, 
where pbo is the central baryon density, Vc the core ra- 
dius, and 13 = 0.4 1.2 (MME99; OM04; Mulchacy ct al. 
2003; Ettori et al. 2004). The baryon profile is also cut 
off at the virial radius. Assuming hydrostatic equilib- 
rium and the NEW halo profile (Navarro, Prenk & White 
1996), Makino, Sasaki & Suto (1998) find that the best- 
fitting /3-model has a core radius about one-fifth of the 
scale radius, rg, of the NEW profile. We choose /? = 0.7 
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Fig. 1. — Halo-modcl mass PS at z = 0. The total PS A'^{k) 
(solid line) is calculated using equation (3) as a weighted sum of the 
dark matter PS A^j(fe) (dashed line), the baryon-dark matter cross 
PS A^^(k) (dotted line), and the baryon PS A.'^^^(k) (dash-dotted 
line). Prom left to right, the thin dash-dotted lines represent one- 
halo contributions to the baryon PS from mass ranges of lO^^M© > 
md > 1O"M0, 1O"M0 > md > W^^Mq, 10^^ Mq > > 
10^"Mq, 10^"Mq > md > lO^M©, and lO^M© > > W^Mq. 



and I'c = 0.22rv to bo the fiducial model and allow the 
parameters to vary. Note that the polytropic gas profile 
(Komatsu & Seljak 2001) is slightly smoother than our 
fiducial model and it gives similar results. 

Fig. 1 shows the fiducial model PS at z = in di- 
mensionless form, i.e. A^(fc) — fc'^P(fc)/27r^. We have 
assumed a flat universe with = 0.3, F = 0.18, 
ug = 0.85, and n = 1, where F is the shape parame- 
ter of the PS, as is the rms value of the density contrast 
within a radius of 8 h~^Mpc, and n is the power spectral 
index. Because, for the fiducial model, the baryon profile 
is smoother than the dark matter profile, the baryon PS 
A^^{k) and baryon-dark matter cross PS A^j(fc) become 
significantly lower than the dark matter PS A^j(fc) on 
small scales {k > 1 h Mpc-^) where the one-halo term 
dominates. At higher redshift, this departure occurs at 
a larger wavenumber. Consequently, the total mass PS 
is lower than what it would be if the hot baryons are 
replaced by dark matter. 

From the break-down of one-halo contributions to the 
baryon PS, also shown in Eig. 1, one sees that each group 
of halos dominates a certain range of scales. The /3-model 
might not apply to halos with md < 10^^ Mq, but since 
they contribute little to the mass PS at fc < 100 h Mpc"-'^ 
- the most important scales to shear PS at I less than a 
few thousand - our conclusions will not be affected. 

3. SHEAR POWER SPECTRUM 

With the Limber approximation and the assumption 
that sources are on a thin slice at z, the shear PS C; in a 
flat universe is given by (Bacon, Rcfrcgier & Ellis 2000; 
Hu 2000; Bartelmann & Schneider 2001) 
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where Hq is the Hubble constant at z = 0, c is the speed 
of Hght, D and Dg are the comoving distance of the lens 
and sources, respectively, a = {1 + z)~^ , and k = l/D. 
We define C; = 2C;/7r, which is the contribution to the 
variance of the deflection angle from logarithmic intervals 
in I. The results are shown in Fig. 2 for the fiducial 
model. Since the shear PS is mostly affected by the linear 
and quasi-linear part of the mass PS, the effect of hot 
baryons on the shear PS is less than a few percent for 
I < 3000. As redshift increases, the same / corresponds 
to smaller wavenumbers, meanwhile the PS is more linear 
at the same wavcnumber. In addition, hot baryons make 
less of a difference in the mass PS at higher redshift. 
Thus, the effect of the ICM on the shear PS reduces 
with increasing redshift. 

The effect of the ICM on the shear PS is large enough 
that it will need to be addressed if the cosmological pa- 
rameter errors forecasted in. e.g.. Song & Knox (2003, 
with Imax = 1000), Hu (2002, with ;,nax = 3000) and 
Refregier et al. (2004, with Zmax = 2 x 10^) are to be re- 
alized. This is clear from the lower panel of Fig. 2 where 
departures outside the shaded region arc larger than the 
statistical error for a fiducial survey with parameter val- 
ues as given in the caption. The statistical error in Ci in 
a band of width Z/4 is given by 

where /sky is the fraction of sky covered, n is the surface 
number density of sources with measurable shapes, and 
7rms is the variance of the "shape noise" . 

From the break-down of the contributions from differ- 
ent mass halos, one sees that the two-halo term in the 
mass PS is the largest contributer at large angular scales, 
and only those halos with rna > lO^^M© are important 
to the shear PS at I around a few thousand. 

Observationally, the (3 parameter and the core radius 
assume a range of values. In fact, the baryon profile 
becomes more compact than the dark matter profile if 
/3 is larger and the core-radius-to-scale-radius ratio is 
smaller. This can drive the baryon-dark matter cross 
PS and baryon PS higher than the dark matter PS, and 
results in an increase in the shear PS. On the other hand, 
the entropy floor (Ponman, Cannon & Navarro 1999) 
can increase the core radius of the baryon profile in low 
mass systems (Tozzi & Norman 2001; Holder & Carl- 
strom 2001) and reduce the total mass PS and shear PS. 

Fig. 3 explores the relative changes in the shear PS 
in the (3-rc/rs space for I = 1000 and 3000. We see 
that given the observational uncertainties about the (3- 
model parameters, hot baryons could pose a significant 
challenge to precision cosmology. 

We are optimistic that the challenge can be met, at 
least for I < 3000, by a combination of hydrodynami- 
cal simulations (which in principle can track hot baryons 
quite accurately), and observations of hot gas via X-ray 
emission and the Sunyaev-Zel'dovich (SZ) effect. For- 
tunately, it is only very massive halos that are impor- 
tant at I < 3000, and these are just the ones easiest 
to study in X-ray and SZ. Planned SZ surveys such as 
South Pole Telescope, Atacama Cosmology Telescope, 
Atacama Pathfinder Experiment and Sunyaev-Zel'dovich 
Array will have sufficient angular resolution and observe 
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Fig. 2. — Upper panel: Shear PS with and without baryons (solid 
lines and circles, respectively) . From bottom to top, the solid lines 
and circles assume sources to be on a thin slice at z = 0.5, 1, and 2. 
The dotted lines represent one-halo contributions to the total shear 
PS Ci a,t z = 0.5. Contributions from halos with mass less than 
are below 10^^^. The fractional difference are shown in 
the lower panel. The dark grey band is the one a statistical error 
on a band of width 1/4 in the sample variance limit for a survey 
with /sijy = 0.25. The light grey band includes shape noise for the 
Zs = 1 PS assuming n = lO/arcmin^ and 7rms = 0.2. 



a sufficiently large number of clusters to guide the devel- 
opment of a statistical model of ICM profiles. Observing 
strategy can also, to some degree, mitigate the difficul- 
ties posed by baryons. Their importance at small scales 
argues for large sky coverage, to reduce sample variance 
(thereby improving the Ci measurement on all scales), 
rather than depth and angular resolution, which reduce 
the shape noise that is only important at small scales. 

4. DISCUSSION AND CONCLUSIONS 

We have extended the halo model to calculate the total 
mass PS of dark matter and hot baryons. Given the 
observed properties of the ICM, we find the total mass PS 
is considerably lower than the dark matter PS on small 
scales. This leads to a few percent reduction of the shear 
PS at 1000 < I < 3000 compared to what it would be 
if baryons traced dark matter exactly. The effect grows 
with increasing I. Further, for I > 3000 an effect (with 
opposite sign, see Fig. 2) also becomes important due to 
cooling/cooled baryons (White 2004). 

So far, we have not considered the ICM (including the 
warm-hot IGM, Dave et al. 2001), which contains roughly 
two thirds of the baryons at ^; = 0. Using hydrodynam- 
ical simulations (TrecSPH, Dave, Dubinski & Hernquist 
1997) that incorporate cooling, heating, star formation, 
and feedback, Zhan (2004) obtains a total baryon PS that 
is lower than what has been shown in Fig. 1. Therefore, 
when all the baryons are accounted for, the effect on 
weak lensing statistics could be even greater. A detailed 
investigation using hydrodynamical simulations will be 
necessary to assess the total baryon effect. In conclu- 
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Fig. 3. — Contours of the fractional difference, Cj/C;, dd — 1) in P^c/ts space. Grey dots in the lower left panel correspond to the observed 
values of /3 and Tc/ts from MME99 and OM04. 



sion, for precision cosmology with weak lensing to live 
up to its promise, we will have to pay attention to the 
modeling of baryons, and not just the dark matter. 
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